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ABSTRACT 

We use the relativistic hydrodynamics code Cosmos++ to model the evolution of the radio 
nebula triggered by the Dec. 27, 2004 giant flare event of soft gamma repeater 1806-20. We 
primarily focus on the rebrightening and centroid motion occurring subsequent to day 20 following 
the flare event. We model this period as a mildly relativistic (7 ~ 1.07 — 1.67) jetted outflow 
expanding into the interstellar medium (ISM). We demonstrate that a jet with total energy 
~ 10^^ ergs confined to a half opening angle ~ 20° fits the key observables of this event, e.g. the 
flux lightcurve, emission map centroid position, and aspect ratio. In particular, we find excellent 
agreement with observations if the rebrightening is due to the jet, moving at 0.5c and inclined 
^ 0° — 40° toward the observer, colliding with a density discontinuity in the ISM at a radius of 
several 10^^ cm. We also find that a jet with a higher velocity, > 0.7c, and larger inclination, 
> 70°, moving into a uniform ISM can fit the observations in general, but tends to miss the 
details of rebrightening. The latter, uniform ISM model predicts an ISM density more than 100 
times lower than that of the former model, and thus suggests an independent test which might 
discriminate between the two. One of the strongest constraints of both models is that the data 
seems to require a non-uniform jet in order to be well fit. 

Subject headings: pulsars: individual (SGR 1806-20) — stars: winds, outflows — hydrodynamics 
— relativity — shock waves 

1. Introduction 

On Dec. 27, 2004 a giant flare (GF) of 7-rays was observed from soft gamma repeater (SGR) 1806-20. 
This GF, the largest 7-ray transient event ever observed from within our galaxy, was comprised of a 0.2 s 
burst so bright that it briefly saturated all 7 and X-ray detectors in the solar system, followed by a several 
100 s fading X-ray transient (Hurley et al. 2005). Such an outburst is thought to derive from a cataclysmic 
reorganization of the surface magnetic field of the magnetar, a neutron star distinguished by its exceptionally 
high surface magnetic field, ^ 10^^ G (Duncan & Thompson 1992), which is the source of the SGR. 

What is more, about one week after the GF a radio nebula was observed at the location of the SGR 
(Gaensler et al. 2005). Over the next couple of months this radio nebula was seen to expand and move across 
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the sky (Taylor et al. 2005; Gelfand et al. 2005). The observations suggest two epochs of evolution. Epoch 
I spanned from the first detection of this radio afterglow about 7 days after the flare event to about 20 days 
post flare. This epoch was distinguished by a steep decay in observed flux (~ t~^''^), a slight motion of the 
flux centroid toward the West, and the detection of polarization of a few percent (Gaensler et al. 2005). Also, 
the nebula was observed to be elongated along the direction of motion with a roughly 2:1 aspect ratio. As 
pointed out by several authors, this epoch had interesting substructure, including a steepening of, or break 
in, the flux decay rate at ~ 9 days as well as variation in amplitude and pitch angle of observed polarization 
(Gelfand et al. 2005; Taylor et al. 2005; Cameron et al. 2005). The observation of polarization is indicative 
of some global order to the magnetic field hosting the synchrotron electrons. 

Epoch II, commencing at about 20 days, was distinct from Epoch I in several ways. First, a rapid 
rebrightening occurred, followed by a shallower flux decay law (~ t~^-^) that persisted for the remainder of 
the reported observations, until ^ 80 days. Secondly the centroid of the flux began to move rapidly toward 
the Northwest. Also, polarization abruptly vanished, with no significant detections during this entire period, 
indicating a separate and distinct population of radiating electrons from that of Epoch I (Gelfand et al. 
2005; Taylor ct al. 2005). 

The detection of proper motion of the emission source strongly suggests some sort of asymmetry. This 
asymmetry could be rooted in i) the ejected material having a particular direction, or it) the surrounding 
environment having a particular shape. While nature might easily incorporate both components, the models 
proposed thus far, including the bulk of our work here, fall into the former category. We have, however, 
briefly explored the latter possibility by simulating a spherical isotropic mass ejection event centered in a 
cylindrical, evacuated cavity. This geometry mocks up the cavity that could be swept out by the SGR's spin- 
down luminosity as it moves through the interstellar medium (ISM). In this model, the source of asymmetry 
would be the unhindered propagation of ejecta down the cylinder, in contrast to the cjccta colliding with its 
wall. Instead, the result we found was that the ejecta progressing down the cylinder, diminishing as is 
quickly dwarfed by the nearly isotropic external shock. Thus any early asymmetry is washed out. As such, 
this model makes the strong prediction that the flux centroid should return to the position of the SGR after 
some time. Since this has not been observed to happen for SGR 1806-20, such a model is disfavored. 

The rest of the paper is organized into four sections. We begin by describing the models we wish to 
investigate and frame the parameter space to be explored with some theory, §2. In particular, we study two 
distinct models for the emission of Epoch II: collisional brightening, §2.1, in which this epoch commences 
when ejecta collides and shocks with a density step in the ISM, and doppler brightening, §2.2, for which 
the ejecta's motion through a uniform ISM, and its progressive shocking and decelaration, is responsible for 
the Epoch II rebrightening. The latter model is most similar to that put forward by Granot et al. (2006), 
except that we assume Epoch I is independent of Epoch II (see §2). We also investigate structured jets §2.3 
which we find to play a surprisingly important and interesting role in fitting the data. We then describe 
our numerical modeling method, §3, including our novel approach to post-processing observations from the 
data, §3.1. We describe our results in §4. While we find that both models studied are viable, the collisional 
brightening model, §4.1, does a better job than the doppler brightening model, §4.3, of fitting the data, 
particularly that of the lightcurve and centroid position. Of particular note, we find that structured jets, 
§4.2, rather than uniform jets, are required to fit the data for both models. Finally, we discuss our results 
and further ways these two models might be tested, §5. 
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2. The Models 

In this paper we model the observations of the radio nebula from the Dec. 27, 2004 event. In particular 
we focus on Epoch II, defined by the rebrightening of radio flux and proper motion of its centroid and assume 
it to be distinct and independent of the emission of Epoch I. We adopt the natural hypothesis that Epoch II 
is the result of the interaction of material ejected during the GF event with an external, interstellar medium 
(ISM) (Gelfand et al. 2005; Taylor et al. 2005; Granot et al. 2006). Our goal is to understand and constrain 
how this interaction might have transpired. In particular, we attempt to address four key issues a) the rapid 
acceleration of the centroid motion, commencing at ~ 20 days, b) the rapid transition, Atpk ~ 10 days, from 
rise to decay of the rebrightening peak, c) the final decay slope oc and d) the observed 2:1 aspect ratio 

of the elongated emission ellipse, with major axis along the direction of centroid proper motion. 

To address these issues we begin by assuming two independent flux components. The first gives rise to 
the Epoch I radio flux, independent of that of Epoch II, and can be described as a stationary point source 
at the location of the SGR, empirically parameterized as a decaying power-law 



Fi{t) = 12 
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-2.7 



Jy , 



(1) 



where the frequency dependence is motivated by Figure 1. We also assume that this flux component is 
present throughout all observations and both epochs, but will be referred to as the Epoch I flux component. 
While invoking equation 1 might gloss over the previously mentioned substructure evident at early times, 
we argue that this approximation is adequate and necessary to address its later interaction with Epoch II. 
For example, while there was a suggestion of proper motion of the radio nebula during Epoch I, it was 
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Fig. 1.— 4.8 GHz data (Gelfand et al. 2005) scaled to 8.5 GHz data (Taylor et al. 2005) by Fa.s/Fs.s = 1.5, 
implying F^, cx This scaled 4.8 GHz is denoted with an asterisk (*) throughout this paper. 
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not monotonic, but instead the centroid returned to effectively zero offset; nor was the motion in the same 
direction as the proper motion of Epoch II, and thus was not obviously directly related. The second flux 
component, designated Epoch II, is entirely derived from the shock interaction of ejecta emitted during the 
GF with the ISM. We wish to test the extent to which this simple two-component model can describe the 
data. 

The assumption of a two-component model is a departure from that put forward by Granot ct al. (2006), 
which proposes that the Epoch I emission component is due to the collision of the ejecta with the ISM at 
some time prior to the first observation at 6.9 days. For the reasons enumerated above, we interpret Epoch 
I to be physically distinct from Epoch II. This assumption is further bolstered by the observation that the 
distance the centroid has moved away from the location of the SGR is correlated with the size of the semi- 
major axis along the same direction, suggesting that the SGR continues to contribute to the flux throughout 
Epoch II. 

Within this two flux component framework, our investigation forks into two branches. The first explores 
collisionally brightened flux in Epoch II, whereby the material ejected in the GF coasts largely unhindered 
until hitting a discontinuous jump or wall in the ISM, which induces a strong shock and rapid flux brightening. 
The second branch is Doppler brightening, whereby the GF ejecta is thought to have intercepted the ISM 
prior to the first observation (~ 7 days) and thus will coast, brightening as it expands r^, until exhausting 
its initial kinetic energy to shocked internal energy, at which time it begins decelerating. This model is most 
like that proposed by Granot et al. (2006) except that they propose the shock due to the initial collision (< 7 
days) might account for Epoch I. Such a prediction does not naturally flow from our simulations; although 
we do model the ejecta colliding with the ISM at an early time, Tyjaii/vo ~ 6 days (Table 1), and self- 
consistently calculate any flux deriving from that collision, we flnd any such flux to be negligible. Efforts to 
increase the ISM wall density, n^uaii, were not successful in generating a non-negligible early flux component, 
but instead would significantly impact the ejecta mass and thus alter the conditions of the later evolution. 
Since the focus of this paper is on the later evolution of Epoch II, we assume the two flux components to be 
independent. 



Now we use the empirical emission model for Epoch I (equation 1) to place constraints on the size and 
shape of the emission of Epoch II. At the time of peak brightness of Epoch II, tpk ~ 30 days, the flux centroid 
was offset from the SGR position by dc ~ 130 mas. This is the mean position of both the flux with zero 
offset, di = 0, from the SGR, Fi{t), and the flux from the shocked, ejected material, Fjj 



From equation 1, Fj{tpk) = 1.2 mjy at 8.5 GHz and using the data by Taylor et al. (2005), the total flux 
was Fj + Fjj = 3.9 mJy. Thus we get an estimate for the observed position of the ejected shock at peak 
brightness (when it is beginning to decelerate) 



2.1. CoUisional Brightening Model 




(2) 
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dc = 1.4 X 130 mas = 0.52 
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Thus we estimate the projected velocity of the ejected material across the sky to be djj/tpk ~ 0.5c. The 
projected distance, d, moved by an object over an observer interval, tobs^ can be related to the material's 
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true initial velocity, vq, and inclination angle, 9i, with respect to image plane {6i > is inclined toward the 
observer) by 

tabs = ( 1 - (wo /c) sin 6*4 ) ^-j- (4) 

Vo COS Oi 

and therefore a relationship between the true velocity and inclination is given by 

Q 

"""^^'^ = sin^, + (c/^;^)cose, 

where v±_ = d/tobs- This equation suggests that if v_\_ = du/tpk ~ 0.5c, then vq ~ v± within 10% for 
0° ^ (^i ^ 60°, but increases rapidly to c for inclination angles approaching 90°. Finally, the characteristic 
true distance from the SGR at which the ejecta peaked in brightness was 

R=^- (6) 

COS^i 

Another interesting constraint is the swiftness with which the Epoch II flux rebrightening increased, 
peaked, and began to decay. This entire transition happened within a time span of Atpk ~ 10 days ~ tpk/3. 
If we assume that the peak of emission corresponds to a self-similar Sedov- Taylor blastwavc reaching its 
deceleration radius, then we generally expect Atpk ~ tpk, i.e. it is difficult to make temporal features shorter 
than the blastwave timescale. In our case we assume the emission derives from a jet with radial opening 
angle, measured from its axis. Thus we may estimate this timescale as simply the timespan between 
when the observer first sees the nearest ("front") edge and the furthest ("back") edge of the jet reach a 
characteristic radius, R. Using equation 4 gives 

— (sin^f — sin^h) = 2— 
c c 

where 9f = 6i + 6o and 6b = di — Oq- Surprisingly, this result is independent of the derived quantities R and 
6i and allows one to solve for the opening angle as a function only of observed quantities 

. f cAtpk \ ono/^0-Sc\ /15kpc\ 

If, as has been suggested (Taylor et al. 2005), the front of the jet moves at twice the velocity of the centroid, 
i.e. v± = 2 X dc/tpk ~ 0.7c, then we expect the opening angle to be 9o ^ 14°. 



Atpk = tb — tf = — {sinOf — sm9b) = 2 — sin^o (7) 



2.2. Doppler Brightening Model 

The rapid peaking of Epoch II can also provide a constraint on the angle of inclination and ejecta 
velocity if we assume that this transition is due only to the ejecta beginning to decelerate as it moves into, 
and sweeps up the interstellar medium. This is distinct from the model described in the previous section 
in which the rapid rise in flux of Epoch II is due to the collision of the ejecta with a discontinuous rise 
in external density. As previously mentioned, due to the self-similar nature of this expansion, we expect 
the time to peak tpk oc rdec/vo ~ Atpk, where the deceleration radius, at which the ejecta has expended 
its expansion energy by shocking the external medium, scales like rdec oc [Eq/ {nextVo)Y^^ , where Eq is the 
ejecta energy and riexi is the density of the external medium. If the observed temporal shortening is assumed 
to be due entirely to Doppler blue shift, then the observation Atpk ~ implies (1 — Vos'mOi) ~ 1/3, 

where tpk = d/ {vq cos 9i). Then using equations 4 and 6 we have 

tobs = il-vo sm9i) « i (9) 
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and assuming the proper motion of the head of the jet is twice the value of the centroid 

r dec cos 9i 2dc , , 

vj_ = = ^ ~ ^-^^ (1^) 

allows us to solve for the inclination and velocity uniquely 

Oi = 70° (11) 
vo = 0.7c . 

This allows us to estimate the deceleration radius, r^ec = 2dc/ cos 9i « 1.7 x 10^'' cm (eqn. 10) at an observed 
distance Zobs = 15 kpc. Using this information and the aforementioned estimate of opening angle, we 
can find the total energy of the cjccta as a function of external, interstellar density; p^xt = nf-xtrripC^ . To do 
so we note that at the deceleration radius the initial kinetic energy, Ek , of the ejecta has been exhausted by 
shocking the external medium 

Ek « -frl,cPe.tc\l - 1)(1 + r(7 - 1)) (12) 

where the shocked medium is boosted by the Lorentz factor 7 = (1 — Vq)"^^"^ and the factor (1 + r(7 — 1)) 
is the specific enthalpy of the shocked external medium (e.g. Anninos et al. 2005) and here we take F = 5/3. 
Since we assume the initial ejecta material is cold, the total energy is related to the mass as Eq = ^Mqc^ 
and the kinetic energy as Ej. = {"f ~ l)MoC^, so Eq — Ek/{1 — I/7) and 

E.»lxlO«(A)'(_Ii^)„,3. (13) 

This is a remarkably simple formula for the total energy of the GF ejecta and is determined entirely from 
kinematic arguments. It is independent of radiative efficiency factors which plague most energy estimations. 
We have numerically simulated a jet with the parameters of eqn. 11, iterating to find the best fit to the 
observations for our adopted radiation parameters (see §3.1.2). Our best fit is J7UHI (Table 1). From the 
table we find 0q = 12° and next = 0.003 cm~^, for which eqn. 13 gives Eq = 0.23 x 10^^ ergs, which is about 
three times lower than the actual value from the table as taken from the simulation. Given the sensitive 
dependence of r^ec on Oi, this difference in total energy measurements corrsponds to ^ 6° uncertainty in 6i. 
This uncertainty is significantly less than the jet opening angle, 60, and thus is within the uncertainty of 
the location of maximum emission, due to Dopplcr effects, on the shock front. Therefore we find the simple 
eqn. 13 to be tolerably accurate and insightful expression for the Doppler brightening model. 



2.3. Structured Jets 

A key point of this study is to explore possible jet structures and their effects on obscrvables. The 
physical motivation for this is clear; nature need not contrive that the material was ejected in the SGR as a 
uniform pancake with a distinct edge. The actual morphology of the ejecta is an open question and, as it 
pertains to gamma-ray burst afterglows, has been studied by several authors (Salmonson & Galama 2002; 
Rossi et al. 2002; Kumar & Granot 2003; Salmonson 2003; Rossi et al. 2004). 

We study three basic jet morphologies, all characterized by an opening angle ^0 measured from the jet 
axis. The "uniform" jet is the simplest, where the ejecta mass density, p, and initial velocity, vq, are constants 
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for the region 6 < do, but are both otherwise zero. The "power-law" jet is constructed by multiplying the 
ejecta mass density by the angle factor 



Ap{e) 



1 



1 - [e/e^f 



(14) 



and the velocity is v{9) = voAp{9), where s is typically between 0.5 — 1.0. The "gaussian" jet multiplies the 
mass density by 



Ag{0) = exp 



^0 



(15) 



and the velocity is v{9) = voAq(9). These three jet types arc indicated by a "U", "P", or "G" respectively 
in Table 1 and the power of the velocity factor, s, is unity unless otherwise given in parenthesis. 



3. Numerical Modeling: Cosmos++ 

In this paper we report numerical simulations using Cosmos++ , a recently completed multidimen- 
sional, radiation-chcmo-magnetohydrodynamic code for both Newtonian and (general) rclativistic flows. 
The capabilities and tests of Cosmos++ are described in (Anninos et al. 2005). This code represents a 
significant advance over our previous numerical code (Cosmos; Anninos & Fragile 2003), with improvements 
in ultra-relativistic shock-capturing capabilities, regularization of the symmetry axis, more accurate energy 
conserving methods, and adaptive mesh refinement (AMR). The mildly relativistic flows considered in this 
paper require accurate relativistic hydrodynamical capabilites and substantial resolution to resolve the nar- 
row, length-contracted shock fronts and flow features. Thus this present work requires of all of these features 
and advances. 

The simulations are performed using special relativistic hydrodynamics and are initialized with a cold 

shell of cjccta mass Mo, and total energy at radius vq from the SGR moving with an initial radial speed 
Vq. Since the ejecta exhibited motion in a particular direction, we take the shell to be an axisymmetric jet 
with a characteristic radial opening angle from the jet axis. This angle, is either a distinct edge to a 
uniform jot. or might be a characteristic angle over which density and/or velocity vary as a function of angle 
from a structured jot axis, as was discussed in detail in §2.3. The shell thickness is constrained to be ^ 1 
light-hour, corresponding to the minimum zone size of the refined grid. If the ejecta is linked to the GF event 
it was likely ejected on a much shorter timescale, perhaps seconds, which is consistent with the idea that, 
since the ejecta appears to have a specific direction, it must have been ejected over a fraction of a rotation 
period, 7.6 s, of the SGR unless it was ejected along the rotation axis. However, there was likely a dispersion 
to the ejecta velocity, or a temperature in the gas, so the shell would likely have expanded significantly over 
the ensuing days as it coasted. Thus our initial shell width is not unreasonable. 

The simulations are initiated on a two-dimensional (r, 9) base grid with a typical resolution: 108 x 24. 
The grid covers the spatial domain 0.01 < r < 0.4 light-years, and < 9 < n/S, w/A or 7r/6, with a reflective 
boundary at ^ = 0. The resolution of the base grid is Ar « 1 light-day. Four to six levels of refinement 
are allowed on this base grid during the evolution. The refinement and derefinement criteria are described 
in Anninos et al. (2005) and use the magnitude, slope, and curvature of the density field to ensure that 
the Lorcntz-contractcd shell of the blast wave is always maximally resolved, while low density background 
regions in the flow are minimally resolved. At maximum refinement, the resolution is Ar « 1 light-hour or 
better. The simulations cover a time span from t = 5 — 10 days to f = 75 — 200 days. 
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3.1. Post Processing 

A key step in assessing the merits of a particular simulation is the extraction of observable attributes 
which can be compared with actual observations. In this case we are interested in simulating observations 
at radio wavelengths. To do this wc create synthetic radio intensity "images" for a specified frequency at 
several observer times. From these can be obtained an observed flux for that frequency and time as well as 
the position and shape of the emitting region; all of these data were observed for the radio afterglow of the 
27 Dec 04 giant flare and can thus be compared. 



3.1.1. Radiative Transfer 

In Cosmos++ the relativistic hydrodynamic equations are integrated along parallel time slices, so that 
for a problem state at a given simulation (laboratory) time, t, all positions, {x,y,z), on the mesh are 
simultaneous in time. In this case the laboratory frame is such that the SGR at the coordinate origin is 
stationary and so, to an excellent approximation, is the terrestrial observer. If, say, the observer is located 
at some large distance Zohs along the z-axis, then he will observe light from the object that is essentially 
parallel to the 2;-axis. Thus the observer's spatial coordinates (aJobs) yobsi -^obs) perceived simultaneously at 
time tohs will transform from those of the simulation simply as 

(iobs, a;obs, yobs, Zohs) = {t- z/c, x, y, z) - (^obs/c, 0, 0, Zohs) + {ro/vo, 0, 0, 0) , (16) 

where the last two terms are arbitrary, constant origin oflBets. The middle term corresponds to look-back 
time from the Earth and will be neglected. The last term accounts for the fact that the simulation is started 
(at t = 0) after the ejecta, moving at some velocity vq, has moved a distance ro from the SGR. 

To construct a flux map at a given observer frequency, v, and time, iobs) we integrate the radiative 
transfer equation through the simulation data, along the z-axis toward the observer 

Iv{tohs, Xohs, 2/obs, ^obs) = j [ji^ - auh)dZohs , (17) 

where is the intensity, ji, the emissivity, and the opacity of the material at a position (tobs, a;obsi Hohs-, Zohs)- 
Finally, the observed flux is given by 

Fv{tohs, Zohs) =11 -^i/(iobs, a;obs, yobs, -^obs) X °bs^ ^"bs ^ ^^g^ 

J J ^ohs 

Distance estimates to SGR 1806-20 have been in the range 12 < Zohs < 15 kpc (Corbel & Eikcnbcrry 2004; 
Figer et al. 2005). See, however, the lower estimate of 9.8 kpc obtained from observations of absorption 
spectra of this OF event (Cameron et al. 2005). We take the high value of this range 

Zohs « 15kpc (19) 

as a worst case with respect to energy requirements, so that our reported flux and energy estimates might 
be relaxed by as much as a factor (9.8/15)^ 0.4, but are not likely to be tightened (see §3.2). 

The emissivity and opacity, which are needed in the observer frame, are transformed from those calcu- 
lated in the proper (rest) frame of the fluid, denoted by a prime, as (see eqns. 4.112 & 4.113 of Rybicki & 
Lightman 1979) 

= a{iy) = — ^ = , (20) 

> ^ j{u)=S'j'{,y') = dy{,^/5) . 
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The proper frequency v' = v/5 transforms with the Doppler factor defined as 

5^[^{l-v.)]-' , (21) 
where is the fluid velocity component along the z-axis (toward the observer), and 7 is the Lorentz factor 

J = [1 - V ■ v\-'/^ . (22) 

Thus wc have a framework to simulate the observed intensity map and flux for any special relativistic 
simulation for which the proper radiation transport quantities a'(;y') and j'{v') can be determined locally 
from the state variables of the simulation. In the next section we describe this for the case of self-absorbed 
synchrotron radiation. 



3.1.2. Self- Absorption Synchrotron Radiation 



The observation of non-thermal (power-law) radio emission from the expanding nebula of this GF sug- 
gests synchrotron radiation from the electrons in the shocked material as the dominant emission mechanism. 
The standard theory of self-absorbed synchrotron emission has been presented in several works (e.g. Rybicki 
& Lightman 1979; Chevalier 1998; Li & Chevalier 1999, and references therein). Our implementation of this 
basically standard model is presented here. The physical picture of this model is that hydrodynamic pro- 
cesses such as shocks will energize the plasma being simulated, thus accelerating its population of electrons 
into a power-law distribution as a function of Lorentz factor 



n(7) = 



(23) 



starting from some minimum Lorentz factor 7„i. A key physical parameter is introduced by assuming that 
only a proportion, Cg, of the total proper energy density, e' = E/^ from the simulation, is represented in the 
electrons. Thus the fraction of energy in electrons is 



•'Tn 



(7 — l)meC?'n(j)d'y . 



Plugging in equation 23, we can solve for the constant C 

v2-P 



c = 



p-2 p-1 



(24) 



(25) 



Another fundamental parameter of the model is the proportion, e^, of the total energy density, e', manifest 
as magnetic field energy density, B^/Stt: 



B = (87rese')^/^ • 



(26) 



The synchrotron emission resTilts from the population of electrons moving in, and being accelerated by 
the iJ-field. The relative orientation, or pitch angle, of an electron's trajectory with respect to the S-fleld 
vector strongly effects the radiation. We treat this pitch angle in an average sense by following the treatment 
of Wijers & Galama (1999). They calculate scaling factors such that the characteristic frequency emitted by 
a distribution of electrons (equation 23) in a field B scales like 



oc XpB 



(27) 
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and the isotropic emission at that characteristic frequency will be 

j'iiy'J a <f>pB . (28) 

We parameterize their functions of electron index p (see Fig. 1 of Wijers & Galama 1999) 

Xp « + 0.267, (29) 

(j)p w 0.53 arctan(p) , 

which give values Xp{P = 2.5) = 0.52 and ^p{p = 2.5) = 0.63. 
So finally, we get for the opacity 

">■) - ^---r-^X^)'"' ^,,BK.,Br"r{^y(^y-"^^y^ . ,30, 

where r() is Euler's gamma function, q the electron charge, and we note the factor (meC^)^'"^ is required to 
correct the dimensional inconsistency of the expression in Rybicki & Lightman (1979). The emmissivity is 
given by 

V^q^C(j)pB fp 19\ /p 1 \ /27rmeCi/'\"^^"^^/^ 



(47r)27rmec2(p+l)^U + 12j^U l^J [ SqXpB ) ' ^^^^ 

The radiative properties of a parcel of shocked material from a simulation depends on the proper 
energy and velocity state variables of the plasma in the simulation, (e', v), as well as the emission modeling 
parameters {p,jrn,^eT('B)- The minimum electron Lorentz factor, 7„i, is typically assumed to be a few. For 
the wavelengths primarily considered here we don't find particular sensitivity to this parameter, and so deem 
it sufficient to set it to 7„i w 2. 

The electron index, p, impacts both the spectral and temporal decay slopes and therefore it is desireable 
to have an independent indication of its value. An intermediate result of this paper is that the 4.8 GHz 
(Gelfand et al. 2005) and 8.5 GHz data (Taylor et al. 2005) are consistent with ^k.g/^s.s ~ 1-5 throughout 
the observational period (Figure 1). Furthermore, the single observation by Taylor ct al. (2005) at 22.5 GHz 
at ~ 23 days is also consistent with the other points when scaled accordingly. This implies F^, oc oc 
j/-(p-i)/2 which, in turn implies 

p « 2.2 , (32) 

which is the value we use throughout this paper. While it is clear that detailed midtiwavclcngth observations 
yield rich spectral substructure (Cameron et al. 2005) it is also the case that the shorter wavelengths, where 
we have the best chance at accurately modeling the emission, can be represented by a single electron index. 



3.2. Comparison with Data 

To fit the simulation to the data one must synchronize the simulation time and observer time by adding 
I'o/va to the simulation time (equation 16) and then determine the normalization, f^r (Table 1), required 
to best match the flux magnitudes. In all of the models presented here and summarized in Table 1 (save 
the defiberate exceptions J5A40SE2, J5A40SE.5, J5A40SEN2) the simulations were caUbrated to match the 
data as nearly as possible, so £n can be thought of as a fine tuning parameter to match the data precisely. 



- 11 - 



Since we find the optical depth of these systems to be negligible, the magnitude of the fliix has the 
dependence 

/ \ / \ (P+l)/4 / \ (l-p)/2 

where e' is the proper internal energy of the fluid in the simulation, Ce and eg parameterize the amount 

of this energy in electrons and magnetic field respectively, and is the final efficiency or normalization 
required to scale the simulated flux to that of the data. As denoted, throughout this paper we assume 10% 
of the energy is expressed in each the electrons and magnetic field. 

A convenient and robust scaling can be applied to the interpretation of these simulation results. As 
discussed in §3, the initial ejecta energy, Eq, is deposited as a cold gas into a fixed volume (presumably 
constrained by the deposition timcscale). Thus Eq cx pvQ. Also, as the cjccta expands and shocks, it will 
heat according to e' oc pv"^. For a class of simulations in which a key hydrodynamic size scale is held fixed, 
i.e. the deceleration radius rdec oc {Eo/n/vQy^^, we have n oc Eq for a given vq and this enforces similar 
hydrodynamic evolution, v ~ vo{rdec/i~)~^^'^ ■ Thus for a given model, the shock heating scales with total 
energy; e' oc Eq. So the simulated fiux scales as 

F. 0. f^S.(i^)(^) (^) (^) E,^^-^y\ (34) 

n oc Eq , 



o.iyVo-iy v8-5GHzy Viskpcy 



where n denotes all external number densities; n^^xt, nint, n^aii- For example, models J5A40S and J5A40SEN2 
differ solely in that their total energy, Eq, and external densities figxti'^intif^waii cire doubled in J5A40SEN2 
from those in J5A40S. Equation 34 then suggests that J5A40SEN2 will have higher flux by a factor 2^-^ w 3.5, 
using equation 32. Indeed, the ratio of factors £n is precisely 1/3.5 so that both simulations match the data 
(see Table 1). In fact, we show below that the resulting flux curves are indistinguishable from one another. 
Not only does this validate the hydrodynamic scaling expected from using a barotropic, P oc , equation of 
state, but also allows for convenient reinterpretation of the data. For instance, the factor, £n, can be used 
to scale Mq, Eg, riext, flint, n^aii by £^^^~^^\ or to scale one of the other terms in equation 34. Similarly, if 
an independent measurement were to constrain the external density or efficiency parameters, one can then 
calculate that effect on the other physical parameters. Alternative distances to the SGR (see discussion at 
equation 19) can be accommodated. The uncertainty in the emission physics parameters can be addressed 
by 

/ , \ -4/(p+5) / X -(p+l)/(p+5) 

^0 OC n OC — — . (35) 



0.1 y VO.l 



3.3. Aspect Ratio 

One metric that can be used to compare the shape of the observed radio nebula to that of simulations 
is the aspect ratio. Observationally, this value is the ratio of minor to major axes of an ellipse fitted to 
the shape of the observed radio nebula. In order to make the most realistic comparison with observations 
possible, we begin with our simulated flux map F^ given by equation 18 and convolve it with a 2D gaussian 
filter, cx exp(— l/2(£i/(j)2) to blur the image (call it F*) where g = +2/^- Herein we set cr = 50 mas 
which is reasonable in that it blurs out small features while not significantly impacting the overall structure 
of the image. The aspect ratio of this blurred flux map, F*, is then the ratio of the root mean square width 
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to height defined by 



where 



and ^ = [x, y] and i = [1, 2] 

4. Results 

We present our results as a series of fits to observations in ascending order of initial ejecta velocity, vq. 
For the collisional brightening models, the velocity range we examine begins at t^o = 0.35c « dc/tpk which is 
roughly the minimum velocity the ejecta could have had and still account for the observed centroid offset, 
dc ~ 130 mas, and continues to vq = 0.5c, for which we find three example solutions. In the case of Doppler 
brightened models, we find higher velocities are required, specifically we consider models with vq = 0.7c and 
0.8c. 



4.1. Collisional Brightening Results 

At the lowest ejection velocity we consider, Vq = 0.35c, we demonstrate the effect on observations of 
inclination angle, 0i, of the jet axis out of the observer view plane. In Figure 2 one can sec the sensitive 
dependence of the Epoch II flux on inclination; the lightcurve is advanced and doppler brightened as 6i is 
increased. Note that all curves fall on the same decay law at late times, as expected due to the waning 
importance of doppler effects as the jet decelerates as noted previously in Salmonson (2003). By adding the 
Epoch I component (equation 1) one gets a total light curve and centroid offset prediction for this simulation, 
as seen at various inclinations in Figure 3. One can see that either the flux or centroid off-set curves can be 
fit reasonably well, but not both. To do so one must go to higher ejecta velocity. 

By increasing the ejecta velocity to vq = 0.377c (7 = 1.08) one can begin to come close to satisfying both 
the flux and centroid offset observations. In Figure 4 are shown three separate models, each individually 
tuned to best match the data. The tuning is rather crude, but a basic trend is highlighted. Increasing 
the density of the wall bounding the ISM, Uwaii, or the region internal to it, riint, causes a more rapid 
brightening of the Epoch II flux bump, which can then be compensated by reducing the inclination, 9i. 
Thus the morphology of the external medium is degenerate with the inclination angle. Note also that model 
J377I25, with a significant internal density, Uint, shows some brightening (and motion of its centroid) prior 
to collision with the wall. This qualitatively matches the data better and thus argues for a non-negligible 
internal medium. 

Next we consider vq = 0.5c. While the vq = 0.377c case nominally satisfied the flux and centroid offset 
constraints, it was a bit shy in satisfying the observed centroid offset at 30 days. Figure 5 shows three 
individually tuned cases which each demonstrate good agreement with the data for vq = 0.5c. Of note is 
the 45° deficit in inclination angle between model J5A355 and the other two. This low viewing angle is 
compensated for by reducing the wall radius, r^uaii and increasing its density n^aii so the jet is seen to 
collide and brighten at the same time as the other models. The other models, J5A40S and J5A40G, contrast 



< y2 > _ < y >2 



JJCF^dxdy 
J J F*dxdy 



(36) 

(37) 
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our two basic structural morphologies of the jet, with angular ejecta mass density variation described by 
equations 14 and 15 respectively. It is clear that each model can be tuned to fit the observations, thus 
making differentiation between the models impossible with present data. Figure 6 shows several different 
variation of model J5A40S, as well as uniform jet model J5A40U. The striking agreement of J5A40SEN2 
and J5A40S speaks to the robust hydrodynamic scaling exhibited in these simulations, as discussed in §3.2. 
Varying the total energy Eq by a factor 0.5 (J5A40SE.5) or 2 (J5A40SE2) gives an idea of how sensitive 
these solutions are to this parameter. The former(latter) case reduces(increases) the deceleration radius and 
thus increases(decreases) relative flux at early time. Finally, J5A40U demonstrates the stark difference in 
the uniform jet evolution to the other, structured jets. This will be discussed further in §4.2. 

In order to compare the observational features of different models, it is useful to construct a synthetic, 
observer-eye-view flux map of the emission, as shown in Figure 7. These images show the position and 
shape of the Epoch II emission with respect to the SGR. Note that all of the non-uniform models exhibit a 
characteristic double brightspot at later time - the leading (top) spot corresponding to the external shock, 
while the trailing spot corresponds to the collision of the flow-focuscd cjccta (§2.3). Figure 8 shows the same 
data as Figure 7, but with the image convolved with a 50 mas standard deviation gaussian point spread 
function. This blurring is meant to simulate the limitation in observable resolution. At this resolution the 
double brightspots are blurred together, but a marginal distinction can be made between the Epoch I and 
II emission regions. 

Taking the images of Figure 8 as being the most realistic, one can compare the predicted aspect ratio 
to that of observations as discussed in §3.3 and shown in Figure 9. The simulations exhibit a signature 
decrease-inflection-increase behavior which is the result of the emergence of the Epoch II flux component 
offset from the SGR (prompting a decrease in aspect ratio) until it surpasses the total flux from Epoch I 
(the inflection) and continues to gain dominance over it (the increase). The simulations show qualitative 
agreement with the data. Indeed, the hint of a decrease-inflection-increase behavior in the data, centered 
at around 30 days or so, is indicative of the veracity of the basic model of the emission being comprised of 
two independent epochs. The rapid decrease in aspect ratio at ~ 30 days corresponds well with the rapid 
brightening of the flux curve and rapid motion of the centroid at this time (e.g. Figure 5), indicating that 
the emergence of the Epoch II flux component was swift, < 10 days. 



4.2. Structured Jet Results 

We find that the uniform jet model fails to produce the t~^-^ decay slope observed for Epoch II, but 
instead decays more steeply. This is unexpected when one considers that the flux scaling for a jet with 
constant opening angle, and expanding with a Sedov- Taylor blastwave profile should go like 

F oc irOa frj' oc r^e'"'' oc r'^-^ oc t-O'^^ , (38) 

where {rOo^r is the volume of radiating electrons with emissivity, j' oc e''^"*'^-'^'* (equations 31 & 32). For a 
strong shock the proper energy goes like e' oc and for a Sedov- Taylor solution, v oc r~^/^, while finally, 
t oc r/v oc r''/"^ (equation 4). Indeed we find these scalings to be good approximations of the simulations - 
except that 6^) is assumed to be constant. In reality, a rarefaction wave will propagate inward (i.e. poleward) 
from edge of the jet at the sound speed, Cg, thus reducing the effective radiating area. This rarefaction wave 
moves like Cg oc \fPjp oc \fd oc v, where P = (F — l)e' and p oc constant for strong shocks. So the jet's 
radiating surface is roughly constant in time as it expands, and so the flux goes like 

F oc rj' oc r-^-^ oc r^ "^^ (39) 
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which is consistent with the minimum (i.e. best fit) asymptotic slopes of uniform jets that we've simulated. 
Therefore, on theoretical and numerical grounds, the uniform jet does not satisfy the late time flux curve. 

By contrast, both the structured and gaussian jets can be made to satisfy the t^^'^ decay slope. The 
reason is a phenomenon we call "flow focusing" , which we find to be quite prevalent in models for which 
the velocity decreases with angle from the jet axis. The basic mechanism is as follows: If the jet is given 
some initial velocity variation v{6) decreasing with angle, the ejecta front will distort as the material in the 
jet core begins to lead that in the wings. Thus, as the ejecta drives a shock into the external medium, a 
reverse shock is also driven back into the ejecta which, due to its distortion, will have obliquity and thus 
impart a poleward velocity into the ejecta. A secondary collision then results as the shocked ejecta collides 
onto the jet axis. This second shock produces a region of high pressure behind the on-axis portion of the 
external, forward shock which, in turn, drives it outward. Thus the on-axis forward shock velocity will be 
sustained for longer than expected in a typical Sedov- Taylor blastwave, and thus can create slower observed 
flux decay rates. An example of the flow-focusing in model J5A355 is shown in Figures 10 and 11. One can 
see a characteristic bowling pin shape to the forward shock at late time due to the added pressure of the 
secondary shock of the ejecta colliding on axis. This simulation can be contrasted with that of uniform jet 
J5A40U in Figure 12, where the lack of obliquity in the reverse shock does not prompt an on-axis secondary 
shock and so the evolution is much more akin to that of a Sedov- Taylor blastwave. The effect of flow-focusing 
on the lightcurve can be seen in Figure 6 where the the uniform jet shows much more rapid flux decay than 
the structured jets. 



4.3. Doppler Brightening Results 

The last set of models we study here do not invoke a discontinuous increase in ISM density with a wall to 
prompt the observed rapid brightening of Epoch II. Instead, these models assiime a sufficiently high Doppler 
factor (equation 21) to brighten and temporally compress the observed flux curve. In general, these models 
require higher velocities and inclination angles in order to achieve sufflcient Doppler factors, as discussed 
in §2. In particular, we take the estimates from equations 8 and 11 and seek a uniform jet model with 
Vo = 0.7c, 6i = 70° and — 12°. The results are shown in Figures 13 and 14. In general, these models 
have a much larger deceleration radius oc {Eo/uexty^^ (Table 1) than the coUisional models. The quahtative 
agreement within errors is acceptable. However, one might ask if it is possible to improve the flt by creating 
a faster rise to peak and a slower decay therefrom. Unfortunately, these two goals arc at odds with each 
other. Obtaining a steeper rise to peak requires a larger vo or 9i in order to increase the Doppler factor and 
expedite the observed evolution. Such an adjustment would, in turn, also increase the decay slope, which is 
already steeper than the observation. 

Figure 15 shows simulated flux maps similar to those of Figure 7. Because the external density, next 
is much lower than that of the previously discussed models (< 1%), the forward shock travels further and 
becomes larger before decelerating. Figures 15 & 13 also demonstrate the level of numerical convergence of 
this model. Models J7UHI and J7UL0 begin identically (requiring high resolution initially to resolve the 
problem features), but the latter progressively sheds its flve layers of mesh reflnement (due to an artificially 
high derefinement threshold in the proper mass density) , having entirely deresolved the forward shock at the 
mid-time of the simulation, until it is entirely on the base mesh at the end of the simulation. The remarkably 
good agreement of the data from these two simulations is indicative of the stability of the results. Figure 16 
compares the proper internal energy e', representative of shock heating, for both problems at the end of the 
run: 0.6 year. The problems exhibit excellent qualitative agreements. 
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In Fig. 17, the images of Fig. 15 for model J7UHI arc convolved with a 50 mas gaussian as was done in 
Figure 8. Once again, at this level of blurring, each epoch's flux component should be resolved. The aspect 
ratio of this model is shown in Figure 18. The more gradual increase in flux from Epoch II in these uniform 
ISM models corresponds to a smoother minimum in aspect ratio (at 20 days) than the sharp inflection 
seen in Figure 9. 

Finally, we explore the possibility of a non-uniform jet plowing into the uniform, low-density ISM to see 
what effect flow-focusing observed in the earlier models might have in this case. The lightcurvc for such a 
model, J8S (Table 1), is shown in Figure 19 for two inclination angles, 60° and 70°, and the aspect ratio for 
these is shown in 18. Because of the relatively long distances traversed by the ejecta, a rather mild velocity 
variation was employed, s = 0.25 (sec equation 14). It is evident that J8S does not show improvement over 
J7UHI in the sharpness of the Epoch II peak, even though a higher initial velocity, vq = 0.8c, was used 
to attempt to increase the Doppler factor. However, it does exhibit a more gradual decay slope, oc t~^-^, 
compared to a t~^''^ for J7UHI, which is consistent with the prediction of equation 39 and is more consistent 
with the observed slope. 

5. Discussion 

In this paper we have numerically modeled the radio nebula of the GF of SGR 1806-20 under the 
assumption that it was comprised of two independent flux components: Epoch I assumed to be centered 
on the SGR and to uniformly decay throughout the observational period; and Epoch II assumed to be due 
to the interaction with the ISM of the mass ejected from the SGR during the GF. In general we find that 
this two component model can satisfy all of the observations considered here: the lightcurve, the centroid 
position, and the nebula aspect ratio. 

We have considered two models for how the Epoch II flux bump evolved. First, the collisional brightening 

model assumes that the rapid brightening of Epoch II is caused by the ejecta colliding with a wall or step 
in ISM density (§2.1). This model is very successful at producing the observed rapid flux brightening which 
precisely coincides with a rapid centroid motion (Fig. 5) and rapid decrease in aspect ratio (Fig. 9). Indeed, 
three distinct models of this class are presented, J5A355, J5A40S and J5A40G, all of which satisfy the data 
despite their varied inclination angles, jet structures, and total energies. Distinct by design, these models 
demonstrate some of the degeneracy in the observable characteristics. However, these models do constrain 
the initial ejecta velocity to be vq ~ 0.5c. They also suggest that if there is a "wall" in the ISM, it must 
be several xlO^® cm from the SGR. If Vq was much below this value, it becomes difficult to satisfy the 
centroid offset data (e.g. Figs. 3 &4), but if vq was much in excess of 0.5c, implying a higher inclination 
angle to keep the centroid data consistent, one finds the flux to brighten more rapidly than is observed and 
a tendency for the late time flux decay rate to become steeper, due to Doppler contraction, than the t~^'^ 
that is observed. Thus we believe these three models to trace out the locus of acceptable parameters for the 
collisional brightening paradigm. 

It is interesting to note that the most successful models in this class had riint differing from riext by only 
a factor of ~ 4. It may well be that models in which riint ~ next can be made to work as well, therefore 
implying that, rather than reaching the terminus of a void, the Epoch II brightening might have occurred 
when the ejecta collided with a small cloud of material. 

The second class of models we consider are Doppler brightened. In these models any c;oIlision with a 
density wall or discontinuity in the ISM is assumed to happen prior to any observations (< 7 days) and 
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so cannot be responsible for the rapid flux brightening of Epoch II. Therefore the brightening is only due 
to expansion of the shock forward of the ejecta as it plows into the ISM. Thus the only way to control 
the shock brightening is by Doppler contraction of the evolution as seen by an observer at high angle of 
inclination (§2.2). Thus these models require higher vq and 6i. The two models presented which best match 
the data arc J7UHI and J8S. Both qualitatively match the data, but miss the finer features in Epoch II: the 
rapid flux brightening coupled with rapid centroid motion and aspect ratio decrease. To improve the flt in 
these respects, the initial ejecta velocity and inclination angle need to be higher than what was used here 
Vq > 0.8c, 6i > 70°. Since increasing these parameters makes the flnal decay slope steeper, a structured jet 
morphology, as in model J8S, becomes more likely. Finally, one could appeal to spectral evolution effects 
such as a hardening spectral index p (equation 32) to foster a more rapid brightening. So although these 
Doppler brightened models do not fit the data as well as the collision-brightened models (with their extra 
degrees of freedom associated with the ISM wall), they might be improved with added physics in future 
work. 

A key difference between these two models is the scale of the external density n,,xt- As can be seen 
in Table 1, the coUisional brightening models J5A355, J5A40S and J5A40G have external densities greater 
than 100 times larger than those of the Doppler brightened models J7UHI and J8S, but all these models 
have a similar total energy, Eq ~ 10^^ ergs. Considering the scaling of eqn. 34, next Eq, an independent 
constraint on the external density might help discriminate between the two models. For instance Granot 
et al. (2006) place a lower bound, next > 10~^, which requires Uext and thus Eq be scaled at least three 
times higher than quoted in Table 1. On the other hand, an upper bound on Ugxt would require next and 
thus Eq to be scaled down for the collisional brightening models, which requires improved shock radiation 
efficiency (eqn. 35) which might become physically untenable. 

Finally, in all of the models we discuss here we find that a non-uniform jet, with a velocity profile 
decreasing with angle from the jet axis, is required to fit the decay slope of Epoch II, t~^'^. As discussed 
in §4.2, this is the result of "flow focusing" behind the external shock, which maintains the strength of this 
shock and thus extends the duration of emission. While some basic tuning was required (see Table 1) to 
fit the data, the general feature of flow focusing seems very generic and ubiquitous for non-uniform velocity 
distributions. To further explore and test this idea, it would be useful to have late time (> 80 days) flux 
measurements to catenate onto those used here and compare with similarly extended simulations. It is very 
reasonable to imagine that jetted outflows generated by nature will be faster at the core than in the wings. 
Thus this effect might appear in other jetted phenomena such as gamma-ray bursts or jetted outflows from 
active galactic nuclei. It will be the subject of future work to explore other instances of this mechanism. 

This work was performed under the auspices of the U.S. Department of Energy by University of California 

Lawrence Livermore National Laboratory under contract W-7405-ENG-48. PCF gratefully acknowledges the 
support of a Faculty R&D grant from the College of Charleston. 
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Fig. 2. — Flux curves of jet collision simulation J35A at four inclination angles, 9i, out of the view plane, 
toward the observer. These are compared with 8.5 GHz data (Taylor et al. 2005) and scaled 4.8 GHz data 
(Gelfand et al. 2005), each with the Epoch I flux component subtracted off. Inclination angle strongly affects 
the timing and brightness of observed flux. The best fit for this jet simulation is 6i = 55°. 
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days 

Fig. 3. — Flux curves (top) and centroid position offsets in milliarcseconds (bottom) for simulation J35A 
seen at four inclination angles, 9i. Compare with Fig. 2 with the same data, but without the Epoch I flux 
component. Jet inclination 9i = 55° most accurately matches the lightcurve. The flux centroid has zero offset 
from the SGR until the shock from the ejecta collision with the external medium brightens. Thus we see a 
rapid acceleration of the centroid position between 20 and 30 days. For inclination 9i = 55°, initial ejecta 
velocity vq — 0.35c is seen not to be sufficient to produce the observed centroid offset, but the qualitative 
motion is evident. 
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Fig. 4. — Flux and centroid offset curves for three simulations, J377I25, J377I30, J377I35. All three match the 
data quite well, but are individually calibrated to match at individual angles of inclination, 6i = 25°, 30°, 35° 
respectively. Model J377I35 is the same, except for '^o, as J35A in Fig. 3. Model J377I30 uses higher densities 
in the wall, riujaii to produce the observed rapid rebrightening at a shallower inclination. Model J377I25 
adds internal density, Uint, thus causing an earlier rise to the Epoch II flux, which better represents both 
the flux and centroid position data. 
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days 

Fig. 5. — Flux and centroid oflFset curves for three simulations at vq = 0.5c: J5A355, J5A40S, J5A40G. These 
simulations are viewed at inclination angles, 9i = —5°, 40° and 40° respectively. Each of these models show 
excellent agreement with the data, despite wide variations in the position of their external wall of material 
and jet structures. These models are typified by ejecta velocity, 0.5c, which allows consistency with centroid 
offset data (bottom plot), a collision with an external wall initiating the rapid flux rebrightening, and a 
structured jet which enables flow focusing and thus reproduction of the late-time t~^^^ flux decay law. 
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Fig. 6. — Flux curves for a range of models for vq = 0.5c. Model J5A40S shows good agreement with the 
data (see Fig. 5). Model J5A40U is a uniform jet calibrated to the same peak flux. The rapid flux decay 
of this model demonstrates the important role of flow focusing of structured jets in producing the observed 
t~^-^ flux decay law at late time. Models J5A40SE2 and J5A40SE.5 compare the effect of doubling and 
halving, respectively, the total ejecta energy. Eg. Finally, for J5A40SEN2, total energy, Eq, and all external 
densities, riext, n-int, n^aii have been doubled. The fact that the flux curve for J5A40SEN2 can be scaled to 
match that of J5A40S demonstrates a robust hydrodynamic scaling. 
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Fig. 7. — Simulated flux maps of four models from Figs. 5 & 6 at two times. For the three structured jets 
one can see the two bright spots characteristic of flow focusing, while the uniform jet, J4A40U, can be seen 
to fade rapidly. Each of the four frames is 0.08 x 0.12 ly^ in size and both halves of each frame are plotted 
on the same color scale. 
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Fig. 8. — Same data as the previous figure 7, but with the Epoch I (eqn. 1) flux component added at the 
position of the SGR (denoted by the red "+") and the image is then convolved with 50 mas gaussian to 
simulate resolution limitations. Thus the total flux from one of these flux maps corresponds to the total flux 
in Fig. 5. One can see the rapidly fading Epoch I SGR flux component and that 50 mas resolution begins 

tn difFfirpntiai.fi thfi f.wn fliiv rnmnnnfint.'j 
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Fig. 9. — Aspect ratios are plotted for the simulated flux maps of the four models shown in Fig. 8. The 
simulated aspect ratio is defined in §3.3. This is plotted against the axial ratio of the ellipse fit by Taylor 
et al. (2005) to the observed emitting region. All simulations exhibit a decrease in aspect ratio as the Epoch 
II flux component increases, and an inflection and subsequent increase when the Epoch II flux exceeds that of 
Epoch I. Model J5A40U shows the most constant aspect ratio after 30 days because, as seen in Fig. 6, Epoch 
II decays rapidly in this model and thus stays comparable in magnitude to Epoch I. It is intriguing that 
the data tends to hint at the same decrease-inflection-increase evolution in aspect ratio as the simulations, 
indicative of two independent epochs. 
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Fig. 10. — Snapshots from simulation J5A355 at four times. On the left of each panel, the proper internal 
energy e' is shown. On the right half log(p) highlights structure and the coordinate velocity field v is 
superposed. The ejecta collides with external medium and its flow is focused toward the jet axis. This, in 
turn, creates a secondary hot-spot on axis that powers the forward shock, creating the characteristic bowling 
pin shape of the forward shock at late times. At time 0.3 year the 15° lineout (see Fig. 11) is shown for 
reference. 
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radius (light— years) 

Fig. 11. — Radial lineouts through simulation J5A355 along two angles of inclination, ^uno, from the polar 
axis. Proper density, energy and velocity are represented for the same four times as shown in Fig. 10. One 
can see the flow focusing of original ejecta material onto the jet axis, ^uno = 0°, at roughly a radius of 0.05 
- 0.06 It-yr, which then drives the center of the jet forward. The forward shock can be seen to progress 
outward more rapidly for 6'iino = 0° along a time sequence of radii (0.045, 0.055, 0.075, 0.09) It-yr, compared 
with (0.041, 0.051, 0.065, 0.075) It-yr for 6'iine = 15°. 
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Fig. 12. — Same plots as in Fig. 10, but for the uniform jet model J5A40U. In contrast to Fig. 10, the uniform 
ejecta is not imparted with a poleward velocity component, and therefore does not exhibit flow focusing but 
instead generates a more spherical external shock. The left half of the final time panel shows the adaptive 
mesh zoning. 
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Fig. 13. — Flux curves of uniform jet expansion simulations J7UHI and J7UL0 as seen from three different 
inclination angles, 0i — 80°, 70°, and 60°. These models only differ in their spacial resolution; J7UHI uses 
5 levels of mesh refinement on top of the base mesh, while J7UL0, after being started at high resolution to 
resolve the initial motion, was deresolved to the base mesh, and thus has zone size 2^ — 32 times larger than 
J7UHI. The fluxes J7UHI were multiplied by 90% (see Table 1) to best match the data, but otherwise the 
agreement of these two resolutions is remarkable. As per the discussion in §2, the case in which vq = 0.7c 
and 9i = 70° basically satisfies the sharp peak of the Epoch II flux bump. However the decay slope of these 
simulations tends to be steeper than the data. The solid black dots designated a, b, c, d are points at 
which flux maps are displayed in Fig. 15. 
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Fig. 14. — Flux and centroid curves for model J7UHI seen at three inclination angles. While the 9i = 70° 
simulation shows general agreement, note that the flux curves show a rounder, broader peak and steeper 
final decay than suggested in the data. The position of the centroid is dominated by the Epoch 11 component 
and is quite sensitive to inclination angle. 
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Fig. 15. — Simulated flux maps of the evolution of the radio nebula at the four times alphabetically indicated 
in Fig. 13. Here the left half of each 0.16 x 0.221y^ panel is from simulation J7UL0, while the right half is 
from J7UHI. The dot at the lower center of each panel indicates the position of the SGR. Both halves use the 
same intensity scale and exhibit close qualitative agreement despite their large diff'erence in resolution. The 
low-resolution (left) simulation shows two sets of discrete concentric bands centered on the jet and the SGR, 
indicating its limitation in spatial and thus temporal resolution, respectively. As shocks bounce between the 
ejected mass shell and the forward shock, their transient flux is evident as broad bands concentric on the 
SGR; a most obvious example of which is seen in panel b). 
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Fig. 16. — A comparison of proper energy densities for simulations J7UHI (right) and J7UL0 (left) at their 
completion time: 0.6 year. They are plotted on the same color scale. One can see the external shock has 
driven upward and expanded from its initial position, ro = 0.01 light-year. Both simulations used the same 
base mesh, shown on the left for reference, but J7UHI allowed 5 levels of refinement on top of this mesh. The 
structure in the flow behind the shock, at radius ~ 1/3 ly, is due to Rayleigh- Taylor instabilities disrupting 
the original disk of ejected mass. 
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Fig. 17. — Model J7UHI for the same sequence of images as in Fig. 15, but with the Epoch I (eqn. 1) flux 
component added at the position of the SGR (denoted by the red "+") and the image is then convolved 
with 50 mas gaussian to simulate resolution limitations. The panel dimensions, 0.16 x 0.251y^, are slightly 
enlarged from those of Fig. 15 to accommodate the blurring. 
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Fig. 18. — Aspect ratios for higher velocity simulations show the same qualitative behavior as those in Fig. 9. 
Because these models exhibit a more gradual brightening of Epoch II (Fig. 13), the aspect ratio decreases 
earlier and has a smoother inflection before increasing. 
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days 

Fig. 19. — Flux and centroid curves for model J8S seen at two inclination angles 60° and 70°. J7UHI is 
shown at 70°, as shown in Fig. 14. Note that, while none of these simulations match the data during the 
flux rise of Epoch II, ^ 20 — 30 days, the structured model J8S gives a final fiux decay slope that is more 
consistent with the observed t~^'^. 
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Table 1. Table of simulations. 



Name 


vq/c 


Bo 


jet (s)(») 










' wall 




'Hnt 


riwall'^^^ 


Figures 




CoUisional Brightening Model 


J35A 


0.35 


24° 


P 


3.56 


3.26 


1.15 


2 


3.15 


1.0 


10-6 




0.6 


2-3 


J377I25 


0.377 


24° 


P 


3.08 


2.77 


1.05 


2 


3.15 


1.0 


0.26 




6.0 


4 


J377I30 


0.377 


24° 


P 


3.08 


2.77 


1.05 


2 


3.15 


1.0 


10-6 




6.0 


4 


J377I35 


0.377 


24° 


P 


3.08 


2.77 


1.15 


2 


3.15 


1.0 


10-6 




0.6 


4 


J5A355 


0.5 


24° 


P 


2.53 


2.35 


0.7 


2 


3.5 


1.04 


0.26 




5.1 


5,7-11 


J5A40S 


0.5 


24° 


P (0.5) 


1.56 


1.5 


0.95 


1 


4.9 


0.53 


0.1 




1.25 


5-9 


J5A40U 


0.5 


12° 


U 


0.2 


0.21 


0.95 


1 


4.9 


0.53 


0.1 




1.25 


6-9,12 


J5A40G 


0.5 


15° 


G (0.75) 


1.09 


1.03 


1.15 


1 


4.9 


0.83 


0.21 




1.3 


5,7-9 


J5A40SE2 


0.5 


24° 


P (0.5) 


3.13 


3.0 


0.4 


1 


4.9 


0.53 


0.1 




1.25 


6 


J5A40SE.5 


0.5 


24° 


P (0.5) 


0.78 


0.75 


2.7 


1 


4.9 


0.53 


0.1 




1.25 


6 


J5A40SEN2 


0.5 


24° 


P (0.5) 


3.13 


3.0 


0.271 


1 


4.9 


1.06 


0.2 




2.5 


6 


Doppler Brightening Model 


J7UHI 


0.7 


12° 


U 


0.6 


0.76 


0.9 


1 


1.1 


0.003 


3 X 10" 


-7 


0.003 


13-19 


J7ULO 


0.7 


12° 


u 


0.6 


0.76 


1.0 


1 


1.1 


0.003 


3 X 10- 


-7 


0.003 


13,15-16 


J8S 


0.8 


20° 


P (0.25) 


1.5 


1.9 


0.8 


1 


1.1 


0.003 


3 X 10" 


-7 


0.003 


18-19 



(")jet type: (U)niform, (P)ower-law, (G)aussian and velocity power, s (=1 unless specified in parenthesis). See §2.3. 

(b) X 10^5 gm 
(°) X 10*6 ergs 

('^^Fine tuning parameter to match simulated flux to data (see §3.2). Multiply Mq, Eo, leit, ^ijntj ^loaH by fj/^'*. 

X IQ-^ light-year 
(■'^^bajryons cm— ^ 



